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INTRODUCTION 


The  extensive  use  of  space  for  platforms  for  communications,  surveillance  such 
as  weather  and  earth  resources,  science  research,  military  objectives,  and  manned  activi¬ 
ties  results  in  a  critical  requirement  for  knowledge  about  the  energetic  particle  environ¬ 
ment  in  space.  Radiation  damage  to  circuits  and  materials,  background  effects  in  sensors, 
hazard  to  personnel,  spurious  effects  in  circuits:  all  are  the  result  of  this  energetic  particle 
environment.  Utilization  of  space  is  continuing  to  increase.  With  this  increase  comes  an 
equivalent  increase  in  the  number  of  personnel  who  have  to  have  basic  knowledge  about 
and  access  to  information  about  the  space  particle  environment.  The  material  here  is  in¬ 
tended  to  provide  an  overview  in  the  areas  of  the  dynamics  of  the  particle  environment, 
trapped  radiation  morphology,  current  trapped  radiation  models,  and  future  modelling 
activity  in  this  field.  The  references  accompanying  this  discussion  can  serve  as  a  conve¬ 
nient  source  for  more  detailed  information  in  this  field. 

A  Gieger-Muller  tube  launched  February  1,  1958  on  Explorer  I  by  Dr.  James  Van 
Allen  is  considered  to  have  produced  the  discovery  of  the  trapped  radiation  belts  sur¬ 
rounding  the  earth,  but  such  a  phenomenon  had  been  predicted  prior  to  the  launch  of  the 
first  artificial  earth  satellite.  A  great  flurry  of  activity  aimed  at  understanding  the  radia¬ 
tion  belts  (sometimes  called  the  Van  Allen  belts )  ensued.  However,  the  myriad  of  mea¬ 
surements  initially  resulted  in  the  acquisition  of  data  rather  than  the  acquisition  of  under¬ 
standing.  Part  of  the  problem  was  the  fact  that  data  were  being  organized  in  terms  of  the 
orbital  parameters  longitude-latitude-altitude,  a  three  parameter  space,  and  the  quantity  of 
data  u'as  insufficient  to  provide  a  sufficient  density  of  data  points  in  any  part  of  space  to 
produce  a  reliable  picture  of  the  particle  population.  In  addition  to  the  three  spatial  pa¬ 
rameters,  there  were  also  energy,  time,  and  species  to  contend  with.  Some  instruments 
made  integral  energy  measurements  (detected  everything  above  a  particular  energy 
threshold,  such  as  gieger  counters  which  counted  any  charged  particle  which  could  pene¬ 
trate  its  window  or  walls),  some  made  differential  energy  measurements,  some  made  uni¬ 
directional  measurements  using  collimators,  others  accepted  particles  coming  from  any¬ 
where  within  a  2k  or  4tc  steradian  angle.  Some  detectors  were  on  spinning  vehicles. 
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other  on  stabilized  platforms.  Since  the  particle  population  is  not  isotropic  (same  inten¬ 
sity  in  all  directions),  further  confusion  was  possible.  The  first  major  step  toward  pro¬ 
ducing  order  out  of  chaos  was  taken  by  Prof.  Carl  Mcllwain  when  he  introduced  a  new 
variable,  L,  based  on  the  second  adiabatic  invariant  of  particle  motion  [1],  which  with  the 
magnetic  field  coordinate  B  formed  a  two-parameter  space  to  replace  the  three  spatial 
coordinates  used  to  describe  satellite  orbits.  The  variable  L  is  discussed  in  the  next  sec¬ 
tion. 

In  1965,  NASA  funded  Dr.  James  Vette  of  The  Aerospace  Corporation  to  produce 
model  electron  and  proton  environments  using  the  data  then  available  from  the  various 
satellites.  All  the  particle  data  sets  available  within  the  US  were  assembled  and  incorpo¬ 
rated  into  model  environments-AEl  (Aerospace  Electron  model  environment  number  1) 
and  API  (Aerospace  Proton  model  environment  number  1)  were  the  result.  These  initial 
models  could  not  be  considered  to  be  more  than  educated  guesses.  Great  effort  was  ex¬ 
pended  in  translating  the  various  data  sets  into  a  common  parameter  space.  (As  noted 
above,  different  variables  controlled  the  various  experimental  measurements  of  the 
trapped  particles.)  Since  with  sparse  data  sets  differential  energy  measurements  could  be 
converted  reliably  into  integral  fluxes  and  unidirectional  measurements  into  omnidirec¬ 
tional  but  not  vice-versa,  the  "lowest  common  denominator”  approach  was  used:  fluxes 
were  specified  as  integral-omnidirectional  fluxes  as  a  function  of  location  in  B,L  space. 
This  is  still  the  mode  used  for  most  trapped  particle  environment  models.  Two  to  three 
orders-of-magnitude  disagreement  between  ostensibly  similar  measurements  were  the 
norm  rather  than  the  exception,  even  in  some  cases  where  the  measurements  were  made 
at  the  same  time  by  different  instruments  on  the  same  satellite.  It  was  realized  that  pro¬ 
ducing  a  reliable  environmental  model  was  going  to  be  a  gargantuan  task.  The  National 
Space  Science  Data  Center  was  set  up  by  NASA  at  the  Goddard  Space  Flight  Center  with 
Dr.  Vette  as  its  director  and  a  major  effort  at  space  environment  modelling  was  begun, 
one  that  continues  at  a  low  level  to  this  day.  The  models  produced  by  those  efforts  will 
be  discussed  later  in  this  presentation. 


TRAPPED  PARTICLE  MORPHOLOGY 


In  this  section,  three  subjects  will  be  discussed:  trapped  particle  motions,  Mcll- 
wain’s  "L"  parameter,  and  the  general  configuration  of  the  radiation  belts.  The  discus¬ 
sion  is  at  an  introductory  level:  no  previous  knowledge  about  the  particle  radiation  envi¬ 
ronment  is  assumed.  The  intent  is  to  provide  sufficient  background  material  about  radia¬ 
tion  belt  morphology  that  the  rest  of  this  tutorial  and  the  other  presentations  in  this  pro¬ 
ceedings  which  deal  with  particle  radiation  effects  will  be  understood  in  context. 

Trapped  Particle  Motions 

An  understanding  of  the  dynamics  of  the  radiation  belts  requires  some  knowledge 
of  the  dynamics  of  an  individual  particle.  The  three  basic  particle  motions  in  the  earth’s 
magnetosphere  with  which  we  are  concerned  are  gyration  or  cyclotron  motion,  bounce, 
and  drift.  The  gyration  is  about  the  local  magnetic  field  line,  the  bounce  motion  is  from 
one  end  of  the  field  line  to  the  other  (one  hemisphere  to  the  other),  and  the  drift  is  around 
the  world  in  longitude.  These  motions  are  a  consequence  of  the  behavior  of  a  charged 
particle  with  forces  acting  on  it  moving  in  a  non-uniform  magnetic  field. 


Figure  1 .  Geomagnetically  Trapped  Particle  Motions 
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In  a  uniform  magnetic  field  B,  a  charged  particle  q  moving  with  velocity  v  expe¬ 
riences  a  force  which  appears  as  an  electric  field  E  at  right  angles  to  both  the  direction  of 
the  field  B  and  the  component  of  the  velocity  vector  perpendicular  to  that  field  (E  =  -q 
VxB  in  vector  notation).  Since  at  each  instant  this  electne  field  tends  to  change  the  direc¬ 
tion  of  the  particle,  which  in  turn  changes  the  direction  of  the  effective  electric  field,  the 
particle  executes  a  circular  path — it  gyrates  about  the  field.  This  coupling  of  forces  also 
results  in  a  complementary  behavior:  a  charged  particle  initially  at  rest  in  a  magnetic 
field  which  has  a  force  imposed  upon  it  will  move  in  a  net  direction  perpendicular  to  both 
the  force  and  the  magnetic  field  (components  of  force  along  the  field  are  ignored  here). 
The  three  basic  motions  of  trapped  particles  in  the  magnetic  field  are  a  consequence  of 
these  force  couplings  and  of  the  fact  that  the  geomagnetic  field  has  curvature  and  inten¬ 
sity  gradients.  Because  of  the  curvature  and  intensity  gradients,  the  particle’s  gyration 
path  does  not  close  upon  itself,  resulting  in  a  drift  motion.  The  motions  are  shown  in 
Figure  1. 

The  direction  of  gyration  follows  the  "right  hand  rule”  for  both  electrons  and  ions; 
since  the  charges  are  opposite  for  the  two  types  of  particles,  the  direction  of  gyration  is 
opposite  for  the  electrons  and  ions,  and  therefore  the  direction  of  drift  is  also  opposite, 
with  electrons  and  negative  ions  drifting  eastward  and  protons  and  other  positive  ions 
drifting  westward.  The  frequency  of  gyration,  called  the  Larmor  frequency,  is  given  by: 

fL  =  Q1/27t  =  -qB/2Ttm0YC  (1) 

where:  q  =  charge  on  the  electron  or  ion 
B  =  local  magnetic  field 
mo  =  rest  mass  of  the  electron  or  ion 
y=  relativistic  mass  ratio  of  the  particle,  1/(1-v2/c2)i/2 
v  =  velocity  of  the  particle 
c  -  velocity  of  light 

Note  that  since  the  frequency  of  gyration  is  proportional  to  B,  it  is  not  constant  along 
the  field  line.  It  is  a  minimum  at  the  magnetic  equator  and  increases  as  the  particle 
moves  away  from  the  equator.  Typical  equatorial  frequencies  at  1000  km  altitude  are 
around  0.5  MHz  for  very  low  energy  electrons  and  about  300  Hz  for  low  energy  protons. 
High  energy  particles  have  lower  gyrofrequencies  because  of  their  greater  (relativistic) 


mass. 


The  First  Adiabatic  Invariant 


If  one  analyzes  the  path  of  a  particle  gyrating  in  a  magnetic  field,  one  observes 
that  it  encloses  a  fixed  amount  of  flux  which  depends  on  the  momentum  of  the  particle 
perpendicular  to  the  magnetic  field.  Under  static  conditions,  the  flux  can’t  "leak"  out  of 
the  path.  This  flux  quantity,  which  is  the  origin  of  the  magnetic  moment  of  the  particle, 
is  invariant  as  long  as  conditions  are  adiabatic;  i.e.,  the  magnetic  field  is  quiescent  and  no 
energy  is  added  to  the  particle.  It  is  called  the  first  adiabatic  invariant.  The  radius  of  gy¬ 
ration  is  inversely  proportional  to  B.  Thus  the  total  flux  enclosed  by  the  path  (flux  den¬ 
sity  times  area)  is  inversely  proportional  to  the  square  of  the  magnetic  field.  But,  if  the 
magnetic  field  increases,  thereby  decreasing  the  radius  of  gyration,  the  perpendicular 
momentum  of  the  particle  has  to  increase  to  conserve  the  magnetic  moment-to  keep  any 
of  the  flux  from  leaking  out  of  the  enclosed  path.  Otherwise  the  flux  area  enclosed  by  the 
gyration  path  would  have  decreased  by  B2  while  the  flux  density  increased  by  only  B. 
This  has  two  interesting  consequences: 

First,  it  produces  the  bounce  motion  of  particles  in  the  geomagnetic  field.  A  par¬ 
ticle  starting  out  at  the  equator  with  a  component  of  velocity  along  the  field  line  will 
travel  a  helical  path  to  lower  altitude  with  the  field  line  as  a  guiding  center.  As  it  moves, 
it  is  moving  in  an  increasing  field.  In  order  to  maintain  a  constant  magnetic  moment,  the 
momentum  component  perpendicular  to  the  magnetic  field  has  to  increase  (in  the  absence 
of  such  an  increase,  the  magnetic  moment  would  be  decreasing  as  1/B).  The  only 
sources  of  energy  to  provide  this  perpendicular  momentum  increase  are  the  magnetic 
field  and  the  particle  velocity.  In  a  quiescent  field,  all  of  the  momentum  increase  is  ob¬ 
tained  from  the  particle  motion  by  converting  momentum  parallel  to  the  field  into  mo¬ 
mentum  perpendicular  to  the  field.  When  this  source  is  exhausted,  the  particle  motion 
parallel  to  the  field  line  is  zero  and  the  particle  is  gyrating  at  a  field  intensity  Bm,  called 
the  mirror  B.  The  gradient  in  the  field  then  reverses  the  process  (called  mirroring)  and 
the  particle  travels  a  helical  path  to  the  other  hemisphere  where  it  again  mirrors  at  a  mag¬ 
netic  field  intensity  Bm  exactly  equal  to  the  previous  one.  The  two  mirror  points  are 
called  conjugate  points  because  they  are  joined  by  the  field  line  which  is  the  guiding 
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center  of  this  helical  motion.  The  particle  actually  spends  most  of  its  time  at  these  mirror 
points. 

The  second  interesting  consequence  which  follows  from  conservation  of  the  first 
adiabatic  invariant  is  acceleration  of  the  particle  by  an  increasing  magnetic  field.  As 
noted  above,  the  only  sources  of  energy  to  maintain  a  constant  magnetic  moment  in  an 
increasing  field  are  the  particle  motion  and  the  field.  If  the  field  is  changing,  the  particle 
may  maintain  its  magnetic  moment  by  increasing  its  perpendicular  momentum  at  the  ex¬ 
pense  of  the  magnetic  field.  Consider  a  particle  that  is  mirroring  so  that  all  of  its  mo¬ 
mentum  is  perpendicular  momentum.  If  the  field  then  increases,  the  field  must  accelerate 
the  particle  to  increase  the  particie’s  perpendicular  momentum  (which  must  be  done  to 
conserve  the  total  flux  enclosed  within  the  particle’s  gyration  path).  Note  that  the  pani¬ 
cle  still  has  the  same  magnetic  moment  after  it  has  been  accelerated  as  it  had  previously, 
but  it  is  now  at  a  higher  energy.  The  reverse  process  also  works:  a  decreasing  field  will 
decelerate  particles.  Geomagnetic  activity  does  both,  producing  radial  displacements  in 
the  process.  If  the  third  adiabatic  invariant,  discussed  below,  is  violated,  the  result  is  ra¬ 
dial  diffusion  (a  process  in  which  particles  initially  on  the  same  field  line  are  transported 
to  higher  and  lower  field  lines)  which  results  in  a  net  increase  of  energetic  particles  in  the 
outer  zone.  These  energetic  particles  then  diffuse  both  inward  and  outward. 

The  Second  Adiabatic  Invariant 

The  bounce  motion  described  above  also  has  an  associated  invariant,  called  the  sec¬ 
ond  adiabatic  invariant.  Basically  it  is  the  total  magnetic  field  energy  contained  within 
the  envelope  of  the  helical  path  between  mirror  points.  (This  function,  like  the  other  in¬ 
variants,  can  be  evaluated  as  a  line  integral.  This  will  be  discussed  in  more  detail  under 
Mcllwain’s  L  Parameter.)  Note  that  if  the  magnetic  field  is  increased,  the  energy  density 
is  increased.  The  radius  of  gyration  is  reduced  to  compensate  (first  invariant,  described 
above),  but  the  path  length  must  also  be  reduced.  Thus  the  mirror  points  must  be  raised. 
Conversely,  if  the  field  weakens,  mirror  points  can  also  be  lowered  and  the  particles 
could  be  lost  into  the  atmosphere.  The  bounce  period  is  only  a  weak  function  of  the 
equatorial  pitch  angle  of  the  particle  (the  angle  between  the  particle  velocity  and  the 


magnetic  field  line  at  the  magnetic  equator)  since  the  particles  spend  most  of  their  time  at 
the  mirror  points.  The  bounce  period  can  be  approximated  by  [2]: 


2n/Q2  =  (4LRe/v)T(y)  (2) 

where:  L  =  Mcllwain’s  parameter 
R*  =  Earth  radius 

y  =  sin  a  where  a  is  the  equatorial  pitch  angle 


and  T(y)  is  given  by  [2]: 


T(y)=1.3802-.31985(y+y,/2) 


The  Third  Adiabatic  Invariant 


Because  the  field  has  a  radial  gradient  and  curvature,  the  radius  of  curvature  of 
the  gyration  or  orbit  about  the  field  line  is  larger  at  the  top  of  the  orbit  than  at  the  bottom 
(top  and  bottom  referenced  with  respect  to  the  earth  radial  direction).  Thus  the  path  does 
not  quite  close  into  a  circle  and  the  next  orbit  starts  slightly  eastward  (for  electrons)  or 
slightly  westward  (for  positive  ions)  from  the  previous  gyration.  This  advance  in  posi¬ 
tion  results  in  a  drift  motion  around  the  earth.  After  one  drift  period  around  the  earth,  the 
particle  will  be  back  at  the  same  location  in  the  field  where  it  started,  provided  the  field  is 
quiescent.  The  locus  of  points  through  which  the  particle  passes  is  called  its  drift  shell. 
The  total  flux  enclosed  by  this  shell  must  again  be  conserved.  It  is  the  conservation  of 
this  flux  function,  or  third  adiabatic  invariant,  that  causes  the  particle  drift  shell  to  close 
after  one  drift  period.  However,  during  the  time  required  for  a  particle  to  drift  around  the 
earth  (Figure  2),  the  magnetic  field  itself  may  change.  During  large  magnetic  storms,  the 
change  can  be  quite  substantial,  up  to  1%.  After  one  drift  period,  the  particle  may  find  it¬ 
self  on  a  different  drift  shell  with  a  different  field  intensity  (and  therefore  a  different  en¬ 
ergy).  The  resulting  violation  of  the  third  adiabatic  invariant  is  the  primary  source  of 
panicle  acceleration  in  the  magnetosphere. 
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The  drift  frequency  is  a  function  of  the  bounce  frequency  (note  the  T(y)  in  Eqn.  2 
shows  up  in  Eqn.  3)  and  can  be  approximately  represented  by  [2]: 

27t/n3  =  -(31V2nY)(Y2-l)(c/Re)2(m0/qB0)D(y)/T(y)  (3) 

where  D(y)  is  given  by: 

D(y)  =  i/i2{5.5208-.4381y-.6397(y  In  y  +  yi/2)} 


and  B0  is  the  value  of  the  earth’s  magnetic  field  at  the  surface  of  the  earth  at  the  equator. 
Figure  2  presents  a  summary  of  approximate  gyration,  bounce,  and  drift  frequencies  for 
electrons  and  protons  as  a  function  of  energy  and  L. 


PROTONS  ELECTRONS 


2  «  8  >  2  «  8 
l  l 


Figure  2.  Gyration,  Bounce  and  Drift  Frequencies  for  Particles  in  the  Earth’s  Magnetic 
Field  [2]. 

Mcllwain’s  L  Parameter 

As  a  result  of  conservation  of  the  first  invariant,  a  particle’s  instantaneous  pitch  an¬ 
gle  as  it  moves  along  a  field  line  can  be  expressed  in  the  form: 
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Sin2  <Xi  /  Bi  =  constant 


(4) 


where:  cx^  is  the  particle  pitch  angle  at  location  i 

Bj  is  the  magnetic  field  intensity  at  the  same  location 


As  a  consequence  of  this  relationship  between  B  and  a,  if  one  knows  the  pitch 
angle  of  a  particle  at  any  point  on  the  field  line  between  the  mirror  points,  one  also  knows 
Bm  (which  is  the  point  at  which  a  =  90°).  The  loci  of  these  mirror  points  as  the  particle 
drifts  around  the  earth  are  two  rings  of  constant  Bm  (one  in  each  hemisphere).  Provided 
one  knows  the  unidirectional  flux  all  along  the  field  line  below  a  point,  the  relationship 
expressed  by  the  equation  above  permits  the  conversion  of  unidirectional  fluxes  to  omni¬ 
directional  fluxes  along  the  same  region  of  the  field  line.  In  a  similar  manner,  one  can 
reconstruct  the  unidirectional  flux  from  the  omnidirectional  flux.  However,  in  practice  it 
is  much  easier  to  measure  the  unidirectional  flux  along  the  field  line  (it  can  be  done  from 
a  single  equatorial  point  by  using  an  instrument  that  scans  in  the  angle  a)  than  to  measure 
the  omnidirectional  flux  distribution  (which  would  have  to  be  done  by  making  measure¬ 
ments  at  closely- spaced  points  all  along  the  field  line). 

As  a  particle  drifts  around  the  earth,  the  conservation  of  the  second  invariant  re¬ 
sults  in  the  particle’s  guiding  center  tracing  out  a  shell  which  connects  the  two  rings  of 
mirror  points.  The  third  invariant  produces  the  result  that  the  shell  is  closed  upon  itself — 
a  particle  remains  on  the  same  shell  as  it  drifts  around  the  earth.  Of  course,  if  the  mag¬ 
netic  field  varies  during  a  drift  period  (or  bounce  or  gyration),  the  adiabat  associated  with 
that  motion  will  no  longer  be  precisely  conserved.  As  previously  stated,  such  violation  of 
conservation  due  to  magnetic  field  fluctuations  results  in  pitch  angle  diffusion,  cross-field 
particle  diffusion,  and  in  changes  in  the  energy  of  the  particles. 

Mapping  of  the  particle  population  in  the  magnetosphere  requires  multi-dimensional 
labelling:  particle  species;  energy;  pitch-angle;  altitude,  latitude,  longitude.  The  task  of 
mapping  the  radiation  environment  is  greatly  simplified  by  reducing  the  three  spatial  co¬ 
ordinates  to  two  magnetic  coordinates,  B  and  L,  which  are  essentially  the  drift  shells  (L) 
and  mirror  rings  (Bm)  described  above.  The  adiabatic  invariant  associated  with  the 
bounce  motion,  I,  is  obtained  by  integrating  the  function  (l-B/Bm)i/2  between  the  mirror 
points.  Since  this  is  awkward  to  do  in  a  nonuniform  field  (the  field  has  to  be  represented 
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by  a  multipole  expansion),  an  approximate  relation  is  derived  which  can  be  related  to  a 
dipole  field:  L3  B  /  M  =  F(i3  B  /  M).  Here  M  is  the  dipole  moment  of  the  geomagnetic 
field.  The  function  F  can  be  calculated  at  a  number  of  points,  generating  an  interpolation 
table.  L  then  becomes  a  simple  calculation  [1],  Note  thiu  for  the  real  magnetic  field,  L  is 
only  an  approximate  representation  of  I,  although  a  sufficiently  accurate  representation 
for  mapping  purposes.  For  a  dipole,  L  =  R,  where  R  is  the  normalized  distance  from  the 
center  of  the  dipole  to  the  equatorial  crossing  of  the  field  line  labeled  "L".  For  our  pur¬ 
poses,  the  dipole  approximation  will  provide  some  understanding  without  belaboring  the 
mathematics: 

R  =  L  cos3  A.  (5) 

B  =  M/R3  (4  -  3  R  /  L)i/2  (6) 

where  R  and  A.  are  the  usual  radial  distance  and  magnetic  latitude  in  a  dipole  field,  M  is 
the  dipole  moment,  and  L  is  Mcllwain’s  parameter.  Note  that  R  and  A.  are  not  sufficient 
to  describe  the  spatial  characteristics  of  the  particle  distributions  since  a  given  particle 
does  not  drift  at  a  constant  R  or  mirror  at  a  constant  X  except  in  a  true  dipole  field  where 
the  azimuthal  symmetry  produces  a  degeneracy.  The  above  expression  shows  that  in  a 
dipole  field,  L  would  correspond  to  the  radial  distance  from  the  center  of  the  earth  to  the 
equatorial  crossing  of  the  magnetic  field  (A,  =  0°). 

The  Radiation  Belts 

The  earth’s  magnetosphere  contains  a  wide  variety  of  charged  particles,  primarily 
electrons  and  protons,  with  energies  ranging  from  the  thermal  (less  than  1  eV)  to  highly 
relativistic  (tens  of  MeV  for  electrons,  BeV  for  protons).  The  ionosphere  contains  a  cold 
plasma,  in  the  1  eV  energy  range,  with  densities  of  the  order  of  lO^/cm3.  The  ionosphere 
is  generally  considered  to  consist  of  the  neutral  and  ion  components  up  to  about  1000  km 
altitude,  with  the  region  above  this  called  the  plasmasphere  (since  the  constituents  there 


are  highly  ionized,  forming  a  plasma).  The  plasmasphere  particle  density  drops  slowly 
until  a  boundary,  called  the  plasmapause,  is  reached  in  which  the  cold  plasma  density 
drops  abruptly  by  about  2  orders  of  magnitude,  from  103  to  10*  per  cm3  to  below  100  per 
cm3.  The  location  of  the  plasmapause  is  local-time  and  magnetic-activity  dependent  but 
is  generally  found  between  L=4  and  L=5  and  generally  follows  the  field  line  to  higher 
latitudes  rather  than  an  altitude  contour.  Beyond  this  region,  hot  plasma  clouds,  with 
temperatures  of  the  order  of  1  to  10  keV  or  more,  are  sometimes  encountered.  The 
plasma  is  heated  by  magnetic  processes  in  the  tail  and  auroral  regions  of  the  magneto¬ 
sphere.  In  the  remainder  of  this  section,  we  will  restrict  our  discussion  to  the  energetic 
particle  populations  >40  keV.  They  will  be  discussed  by  location  or  zone  and  by  species. 

When  Prof.  Van  Allen’s  detector  passed  through  the  magnetosphere  in  a  radial  di¬ 
rection,  the  count  rate  in  the  detector  first  increased,  then  decreased,  and  then  increased 
again.  Two  distinct  zones  of  trapped  radiation  were  being  traversed.  These  were  named 
the  Van  Allen  belts  and  are  commonly  referred  to  as  the  inner  zone  and  the  outer  zone , 
with  a  region  known  as  the  slot  separating  them.  Only  for  electrons  are  these  zones  dis¬ 
tinct.  Figure  3,  which  presents  data  from  a  period  when  copious  fluxes  of  fission  elec¬ 
trons  were  still  present  from  the  Starfish  nuclear  explosion  in  space,  shows  the  minimum 
between  the  inner  zone  and  outer  zone.  The  inner  zone,  which  is  generally  considered  to 
cover  the  region  1.0  <  L  <  2,  has  a  peak  in  flux  intensity  at  about  L=1.5  for  1  MeV  elec¬ 
trons.  The  region  2.0  <  L  <  2.8  is  generally  considered  the  slot  region  where  magneto- 
spheric  processes  result  in  a  low  intensity  of  electrons  during  magnetically-quiet  periods. 
The  process  which  removes  the  electrons  is  a  resonant  interaction  between  energetic 
electrons  and  whistler-mode  (right  circularly  polarized)  electromagnetic  waves.  The  in¬ 
teraction  results  in  some  of  the  particles  being  scattered  to  lower  angles  (relative  to  the 
magnetic  field  line)  such  that  their  new  mirror  points  are  within  the  atmosphere.  The  at¬ 
mosphere  absorbs  them.  Thus  the  slot  region  normally  contains  relatively  low  fluxes  of 
particles.  At  times  of  large  magnetic  storms,  the  slot  can  be  refilled  and  quite  high  flux 
levels  can  be  observed  there  for  a  few  days.  The  location  of  the  slot  is  quite  variable: 
during  large  geomagnetic  storms,  the  minimum  between  the  inner  and  outer  zone  can  be 
very  narrow  and  may  be  displaced  to  a  low  L  value,  even  centered  as  low  as  L=2.0.  Im¬ 
mediately  after  a  storm,  the  slot  may  be  completely  filled  with  electrons  and  so  does  not 
exist.  An  extensive  discussion  of  these  dynamics  is  available  elsewhere  [3],  The  outer 
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electron  zone,  which  is  extremely  variable,  typically  peaks  around  3.5  <  L  <  4.0.  In  the 
outer  zone,  significant  fluxes  of  electrons  with  energies  in  excess  of  5  MeV  are  observed 
after  major  magnetic  storms.  During  extended  quiet  periods,  the  outer  zone  may  almost 
disappear  at  high  energies.  The  difference  in  flux  intensity  from  minimum  to  maximum 
may  be  as  high  as  5  orders  of  magnitude.  See,  for  example,  Figure  2  of  [3]  or  Figure  7  of 
[4]. 


Figure  3.  Inner  and  Outer  Zone  Electron  Belts.  The  numbers  on  the  contours  are  the 
logio  of  the  integral  omnidirectional  flux,  L*.  The  intense  inner  zone 
was  due  to  the  Starfish  fission  source.  For  times  subsequent  to  1968,  the 
maximum  inner  zone  intensity  is  several  orders  of  magnitude  lower. 

The  proton  environment  is  sometimes  separated  into  two  constituents  also,  but  in 
this  case  the  separation  is  done  on  the  basis  of  the  energy  of  the  protons.  The  same  re¬ 
gion  of  space  that  constitutes  the  inner  zone  for  electrons  contains  very  energetic  protons, 
some  with  energies  in  excess  of  200  MeV,  which  are  also  the  result  of  the  cosmic-ray 
albedo  neutron  decay  described  in  the  section  on  inner  zone  electrons. 

There  is  another  source  of  trapped  protons.  Low  energy  protons,  from  either  a 
solar  wind  or  ionosphere  source,  are  accelerated  similarly  to  the  energetic  electrons. 
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Some  of  the  initial  acceleration  for  ionospheric  electrons  is  produced  by  electric  fields, 
especially  in  the  auroral  zone.  The  low  energy  proton  belt  can  be  considered  to  be  com¬ 
posed  primarily  of  protons  with  energies  below  10  MeV.  These  lower-energy  protons  are 
present  in  both  the  inner  and  outer  zones.  There  is  no  slot  such  as  occurs  for  the  elec¬ 
trons.  The  peak  in  proton  flux  intensity  depends  on  the  energy,  with  lower  energy  pro¬ 
tons  peaking  farther  out.  The  1  MeV  proton  flux  is  at  its  maximum  at  about  L  =  3. 

TRAPPED  PARTICLE  POPULATIONS 

In  this  section,  we  will  briefly  describe  the  four  major  components  of  the  magne- 
tospherically-trapped  energetic  particle  population— the  inner  and  outer  electron  zones, 
the  energetic  proton  belt,  and  the  low  energy  proton  population.  The  sources  of  the  parti¬ 
cles  and  their  general  flux  levels  as  a  function  of  L  and  energy  will  be  discussed.  Solar 
flare  protons,  which  are  discussed  in  the  later  section  on  modelling,  will  not  be  discussed 
in  this  section  because  the  contribution  of  flare  protons  which  become  trapped  to  the  flu- 
ences  observed  in  most  orbits  is  negligible.  Rare  protons  in  the  polar  region  may  be  a 
major  concern  for  some  satellites,  but  those  particles  are  not  trapped. 

Inner  Zone  Electrons 

The  source  of  the  inner  zone  electrons  is  a  combination  of  cosmic-ray  albedo 
neutron  decay  (CRAND)  and  radial  diffusion  through  the  slot  from  the  outer  zone.  In  the 
CRAND  mechanism,  cosmic  rays  interact  with  air  molecules  in  the  upper  atmosphere, 
producing  energetic  neutrons,  some  of  which  escape  back  into  space.  Since  neutrons  are 
uncharged,  they  cross  magnetic  field  lines  unimpeded.  However,  some  decay  while  still 
in  the  magnetosphere  and  the  decay  products,  an  electron  and  a  proton,  are  charged  and 
so  become  trapped.  The  end-point  energy  of  the  electron  in  neutron  decay  is  slightly  un¬ 
der  1  MeV.  The  contribution  of  the  neutron’s  velocity  to  the  electron’s  energy  is  small. 
As  a  result,  there  are  few  electrons  with  energies  in  excess  of  1  MeV  in  the  inner  zone. 
Electrons  with  higher  energies  are  present  in  small  numbers,  especially  above  about 
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L=1.65  after  large  magnetic  storms,  but  can  be  ignored  as  a  hazard  to  space  systems  ex¬ 
cept  for  their  background  effects  in  sensors. 


Inner  zone  electrons  below  about  1000  km  ha  -/e  lifetimes  that  are  primarily  de¬ 
termined  by  the  scale  height  of  the  atmosphere.  During  solar-active  periods,  the  in¬ 
creased  scale  height  results  in  a  reduced  lifetime  and  lower  average  fluxes.  This  is  re¬ 
flected  in  the  models  by  having  a  solar  maximum  and  a  solar  minimum  version.  Farther 
out  in  the  inner  zone,  electrons  are  quite  stable,  with  typical  lifetimes  of  400  days  [5]. 
Principal  loss  mechanisms  are  probably  any  or  all  of  the  following:  radial  diffusion  into 
the  atmosphere  (violation  of  the  second  and  third  invariants  caused  by  magnetic  storms); 
interaction  with  whistler-mode  waves  produced  by  lightning  strokes  (the  resonant  inter¬ 
action  between  these  waves,  also  known  as  cyclotron  waves,  and  the  electrons  results  in  a 
lowering  of  the  electron  pitch  angle,  causing  it  to  be  absorbed  by  the  atmosphere  at  the 
end  of  the  field  line);  interaction  with  VLF  waves  from  ground-based  transmitters. 

The  order-of-magnitude  of  the  electron  fluxes  at  L=1.45  in  the  inner  zone  is  as 
follows:  >  108  for  Ee  >  0.1  MeV;  >  10<s  for  Ee  >  1  MeV;  >  105  for  Ec  >  2  MeV.  The 
numbers  represent  the  integral,  omnidirectional  fluxes  cm-2-sec-i.  Below  about  L  =  1.55, 
the  fluxes  are  quite  stable,  with  little  variation  being  observed  over  the  solar  cycle  [6]  ex¬ 
cept  for  altitudes  below  1000  km  where  atmospheric  effects  are  observed.  Above  L=  1.6, 
major  magnetic  storms  inject  electrons  with  energies  up  to  at  least  1.2  MeV  [3].  Figure  4 
shows  the  equatorial  omnidirectional  inner  zone  electron  flux  intensities  as  a  function  of 
L  and  energy. 


Outer  Zone  Electrons 


The  outer  zone  electrons  originate  either  as  solar  wind  electrons  in  the  tail  of  the 
magnetosphere  or  as  ionospheric  electrons  at  high  latitudes  which  are  accelerated  up  the 
field  lines.  Magnetic  field  fluctuations  cause  them  to  be  diffused  radially  inward,  ener¬ 
gizing  them.  The  acceleration  is  a  consequence  of  the  conservation  of  the  first  adiabatic 
invariant  coupled  with  violation  of  the  third  invariant,  discussed  earlier,  by  magnetic  ac¬ 
tivity.  As  the  particles  are  transported  to  field  lines  deeper  in  the  magnetosphere,  the  in- 
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crease  in  the  average  field  intensity  has  to  be  compensated  by  an  equivalent  increase  in 
particle  momentum,  or  energy.  The  various  fluxes  peak  at  different  locations  in  the  outer 
zone  for  different  energies,  with  the  higher  energies  peaking  at  lower  L.  Representative 
outer  zone  fluxes  are  of  the  order  of:  >  10*  for  Ee  >  C.l  MeV,  L  =  6;  >107  for  Ee  >  1 
MeV,  L  =  5;  >  105  for  Ee  >  4  MeV,  L=4.  Units  are  as  described  for  the  inner  zone 
fluxes. 
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Figure  4.  Inner  Zone  Electron  Figure  5.  Outer  Zone  Electron 

Flux  Intensities  [6]  Flux  Intensities  [7] 

Outer  zone  electron  fluxes  are  highly  variable,  with  increases  at  a  given  energy  on 
a  given  L  shell  being  as  great  as  5  orders  of  magnitude  in  less  than  a  day.  These  large  in¬ 
creases  are  caused  by  major  magnetic  storms,  where  D„  >  -150y.  D„  is  a  global  magnetic 
field  disturbance  index  which  is  generally  responsive  to  low  latitude  variations  caused  by 
a  magnetospheric  ring  current.  This  ring  current  is  composed  of  low  energy  ions  acceler- 
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ated  by  the  magnetic  storm.  Typical  decay  constants  for  outer  zone  electrons  are  of  the 
order  of  10  days.  In  addition  to  the  radial  diffusion  of  particles  caused  by  magnetic 
storms,  they  also  cause  pitch-angle  scattering  of  the  particles.  Thus  panicles  which  were 
previously  stably  trapped  on  a  field  line  (had  mirror  points  that  were  above  the  atmo¬ 
sphere  everywhere  along  their  drift  paths),  can  be  perturbed  so  that  they  now  mirror 
within  the  residual  atmosphere  below  100  km  at  some  point  along  their  drift  path.  At  this 
altitude,  the  atmosphere  absorbs  the  panicles.  A  low  altitude  satellite  which  is  normally 
below  the  trapped  radiation  zones  (except  when  traversing  the  South  Atlantic  Anomaly) 
may  suddenly  find  itself  bathed  in  large  fluxes  of  energetic  electrons  at  midlatitudes 
when  it  encounters  these  particles  which  show  up  low  on  the  outer  zone  field  lines 
(sometimes  called  the  horns  of  the  outer  zone).  The  South  Atlantic  Anomaly  is  a  region 
of  anomalously  low  magnetic  field  strength.  Since  particles  mirror  at  a  constant  Bm,  they 
attain  their  lowest  mirror  altitude  in  the  South  Atlantic  Anomaly. 


Energetic  Protons 


The  source  of  the  energetic  protons  which  are  present  in  the  inner  zone  is 
CRAND,  the  mechanism  mentioned  previously.  They  are  quite  stable,  with  minor  varia¬ 
tions  in  intensity  occurring  at  low  altitudes  due  to  variations  in  the  atmospheric  density 
due  to  solar  activity.  Typical  intensities  are  of  the  order  of:  >  104  for  Ep  >  100  MeV  and 
>  103  for  Ep  >  300  MeV,  both  at  L  =  1.45.  Due  to  the  secular  variation  in  the  magnetic 
field,  a  very  slow  decrease  which  may  be  an  indication  that  the  earth’s  field  will  undergo 
a  reversal  in  the  geologically-speaking  near  future  (104  years?),  the  energetic  proton  envi¬ 
ronment  is  also  exhibiting  a  small  decrease  (the  decreasing  field  intensity  is  driving  the 
protons  into  the  atmosphere,  again  due  to  conservation  of  the  adiabatic  invariants).  A  se¬ 
rious  problem  in  particle  modelling  due  to  this  secular  decrease  in  field  intensity  will  be 
discussed  in  the  section  on  models. 
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Low  Energy  Protons 


The  low  energy  protons  with  which  we  are  concerned  here  are  in  the  0.5  to  5 
MeV  range,  since  there  are  large  fluxes  of  such  particles  in  both  the  inner  and  outer  zone 
and  they  have  significant  materials  effects.  Particles  with  these  energies  can  originate  in 
a  number  of  sources:  radial  diffusion  and  energization  of  solar  wind  particles  which  en¬ 
ter  the  geomagnetic  tail,  similar  to  the  outer  zone  electrons;  ionospheric  acceleration  up 
field 


Figure  6.  High  Energy  Protons  in  the  Inner  Zone  [8] 


lines,  with  subsequent  radial  diffusion  and  acceleration;  direct  access  of  solar  flare  pro¬ 
tons.  Typical  intensities  in  the  outer  zone  are:  >  108  for  Ep  >  0.1  MeV;  >  10?  for  Ep  >  1 
MeV;  >  105  for  Ep  >  10  MeV;  <  102  for  Ep  >  100  MeV.  Again,  these  are  omnidirec¬ 
tional,  integral  fluxes  in  units  of  cm-2-sec-i.  While  the  fluxes  are  subject  to  variation  due 
to  magnetic  storm  activity,  the  variations  are  much  smaller  than  for  electrons.  The  pri¬ 
mary  loss  mechanisms  are  deenergization  through  collisions  with  the  residual  atmosphere 


21 


and  charge-exchange,  which  results  in  an  energetic  neutral  particle  which  is  not  trapped 
by  the  magnetic  field. 


CURRENT  STATUS  OF  PARTICLE  MODELS 


In  this  section,  we  will  cover  the  currently  recommended  particle  models,  their 
ranges  and  estimates  of  accuracy,  discuss  briefly  their  sources,  availability,  future  mod¬ 
elling  plans,  and  requirements  for  additional  data.  A  more  extensive  discussion  of  the 
modelling  efforts  at  the  National  Space  Science  Data  Center  is  presented  elsewhere  [9]. 


Inner  Zone  Electron  Models 


The  current  NSSDC  models  which  provide  useful  inner  zone  (L  =  1.2  to  2.4) 
electron  data  are  AE5  [6],  for  solar  minimum,  AE6  [10],  for  solar  maximum,  and  AE8 
[9]  for  either  solar  minimum  or  maximum.  The  energy  range  of  these  models  is  from 
0.04  MeV  to  5  MeV,  although  present  techniques  can  not  make  reliable  measurements  of 
electrons  with  energies  above  2  MeV  below  about  L=1.55  in  the  inner  zone.  They  are 
empirical  models,  being  based  on  in-situ  measurements  of  the  fluxes.  To  produce  the 
models,  data  are  acquired  from  investigators,  corrected,  edited,  averaged,  interpolated, 
extrapolated,  etc;  in  other  words,  it  is  all  very  thoroughly  massaged.  Some  estimate  of  its 
reliability  is  also  made,  but  such  a  step  is  very  subjective.  In  general,  the  temporal  cov¬ 
erage  for  any  data  set  is  six  months  or  longer.  The  measurements  were  made  between 
1964  and  1977  (although  data  up  to  1968  include  Starfish  electrons).  Starfish  contami¬ 
nation  has  been  removed  from  the  data  in  generating  the  above  models.  Newer  data  are 
available  for  incorporation  into  models  [9].  The  accuracy  of  the  models  is  very  good, 
better  than  a  factor  of  two,  for  energies  below  1  MeV  and  L  <  1.65.  Above  L=1.65,  the 
variability  of  the  flux  levels  themselves  produce  uncertainty.  Above  2  MeV,  the  fluxes 
are  extrapolations  and  their  accuracy  is  unknown. 

The  models,  and  codes  for  running  them,  are  available  from  NSSDC.  In  addition 
to  a  tape  format,  the  codes  are  also  available  as  files  on  a  VAX  at  NSSDC.  The  VAX  is  a 
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SPAN  (Space  Physics  Analysis  Network)  node.  Thus,  anyone  with  access  to  SPAN  or 
other  networks  which  can  connect  to  SPAN  (such  as  TELENET,  ARPANET,  BITNET, 
etc.)  can  access  these  files  to  download  them  to  their  home  computer  via  the  network. 
Dr.  James  Green  is  in  charge  of  SPAN  at  NSSDC.  Alternatively,  the  codes  can  also  be 
run  interactively  on  the  NSSDC  VAX  at  no  charge  to  the  user  (other  than  his  own  costs 
of  connecting  to  the  SPAN  network).  An  alternate  interactive  resource  is  EnviroNET 
[11],  which  is  resident  on  a  MicroVAX  at  GSFC  and  is  accessible  as  the  SPAN  node 
ENVNET.  To  access  EnviroNET,  the  user  name  ENVIRONET  and  password  HEN- 
NIKER  are  used. 

Future  plans  for  the  inner  zone  modelling  activity  are  to  incorporate  newer  data 
bases  and  also  perhaps  include  some  storm-time  dynamics.  It  may  not  be  possible  to  ac¬ 
complish  the  latter  task  with  presently  available  data  bases,  but  the  CRRES  mission  [12] 
has  as  one  of  its  objectives  the  acquisition  of  the  data  required  for  producing  dynamic 
particle  models. 


Outer  Zone  Electron  Models 


In-situ  electron  flux  data  in  the  outer  zone  are  far  from  satisfactory  for  generating 
electron  models.  Most  of  the  data  used  are  extrapolations  in  both  energy  and  altitude.  At 
the  geosynchronous  region,  measurements  up  to  1.7  MeV  have  been  available  from  ATS- 
1.  Near  the  equator,  S3  provided  measurements  from  about  L=2.5  to  L=6,  but  only  up  to 
300  keV.  OGO-5  had  an  electron  channel  at  2.7  MeV  but  the  satellite  orbit  inclination 
was  27*  so  it  made  no  measurements  near  the  equatorial  region.  All  other  satellites 
which  traversed  the  equatorial  region  at  high  altitude  either  had  no  high  energy  electron 
measurements  (Ee  >  1.5  MeV)  or  the  energy  threshold  and  detector  efficiency  were  not 
known  with  sufficient  accuracy  to  be  usable  in  modelling.  The  source  of  the  orders-of- 
magnitude  discrepancies  seen  in  comparisons  of  energetic  outer  zone  electrons  (e.g.,  Fig¬ 
ure  12  of  [4])  are  these  uncertainties.  All  other  sources  of  data  used  in  the  models  are 
extrapolations  of  measurements  made  low  on  the  field  line. 

The  current  NSSDC  outer  zone  models  (L  >  2.4)  which  provide  useful  results  are 
the  following:  AE7-Lo  and  AE8  for  solar  minimum  and  AE7-Hi  for  solar  maximum  or 
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long  duration  missions  (>  5  years).  For  geosynchronous  satellites,  another  model  is  still 
relatively  valid— AE3  [13].  For  long  term  missions,  AE7-Hi  is  probably  accurate  to 
within  a  factor  of  two,  especially  for  L  <  6  and  E  <  5  MeV.  However,  if  the  mission  in¬ 
cludes  the  period  a  year  or  two  following  the  sunspot  maximum  when  the  magnetic  storm 
activity  is  greatest,  AE7-Hi  will  err  on  the  low  side  (actual  integrated  fluences  can  be  ex¬ 
pected  to  be  greater  than  the  model  prediction).  AE7-Hi  was  generated  in  response  to 
criticisms  that  the  earlier  models,  AE4  and  AE6,  were  deficient  in  high  energy  electrons, 
and  seriously  so.  In  fact,  in  those  models  the  energy  spectra  were  truncated  at  5  MeV,  as 
is  also  done  in  AE8.  A  comparison  of  the  models  with  in-situ  data  [4]  showed  that  the 
models  predicted  fluxes  that  were  low  by  about  a  factor  of  three,  but  almost  the  entire  de¬ 
ficiency  was  in  electrons  >1.5  MeV.  The  result  was  a  prediction  of  dose  in  heavily 
shielded  satellite  components  that  was  low  by  an  order  of  magnitude.  AE7-Hi,  which 
truncates  the  energy  spectrum  at  7.5  MeV,  has  also  been  criticized  for  truncating  the 
spectrum,  since  electrons  with  energies  up  to  10  MeV  have  been  measured  at  geosyn¬ 
chronous  orbit  [14].  Figure  7  shows  the  equatorial  flux  contours  as  a  function  of  L  and 
energy  which  are  contained  in  AE7-Hi  [15]. 


Figure  7.  AE7-Hi  Equatorial  Electron  Flux  Contours  [15]. 

AE8  exists  in  two  forms,  AE8MIN  and  AE8MAX,  which  are  supposed  to  repre¬ 
sent  the  environment  during  solar  minimum  and  maximum.  However,  both  are  truncated 
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at  5  MeV  and  cannot  properly  mode)  the  solar  maximum  period  when  large  fluxes  of 
very  energetic  electrons  appear.  AE8  uses  a  single  outer  zone  model  and  uses  AE5  and 
AE6  solar  minimum  and  solar  maximum  inner  zone  models.  The  major  difference  be¬ 
tween  AES  and  AE7-Hi  is  in  the  high  energy  electron  dux  at  around  L  =  4.  The  AE7-Hi 
modei  contains  about  a  factor  of  two  greater  flux  at  3  MeV  and  about  a  factor  of  10 
greater  flux  at  4.5  MeV.  At  solar  minimum,  there  are  relatively  few  energetic  electrons 
and  any  of  the  later  models,  AE7-Hi,  AE7-Lo,  AE8MIN,  and  AE8MAX  are  satisfactory 
for  any  use  except  calculating  background  rates  in  detectors.  The  models  are  probably 
accurate  to  a  factor  of  three  for  dose  calculations. 

As  in  the  inner  zone  case,  future  plans  are  to  incorporate  new  data  in  an  attempt  to 
increase  the  accuracy  of  the  models  and  to  include  a  dynamic  model  (addressing  both  the 
prediction  of  electron  increases  due  to  magnetic  storms  and  the  evolution  of  the  fluxes 
after  the  increases).  The  Radsat  portion  of  the  CRRES  mission  [12]  has  two  outer  zone 
electron  modelling  goals:  obtain  data  in  the  outer  zone  nearer  to  the  equator  than  has 
been  done  previously  in  order  to  reduce  the  amount  of  extrapolation  that  has  been  neces¬ 
sary  for  model  generation;  obtain  a  data  base  which  can  be  used  to  generate  a  dynamic 
model. 

Proton  Models 

The  current  proton  models  are  AP8MIN  and  AP8MAX,  again  representing  the 
solar  minimum  and  maximum  periods.  The  major  effect  of  the  solar  variation  is  the  vari¬ 
ation  in  atmospheric  density  at  lower  altitudes:  at  solar  maximum,  the  higher  scale  height 
of  the  atmosphere  decreases  the  energetic  proton  fluxes.  The  models  are  probably  accu¬ 
rate  to  50%  or  better.  They  cover  the  energy  range  from  100  keV  to  400  MeV  and  the  L 
range  from  1.17  to  7.  The  data  were  obtained  during  the  same  time  period  that  the  inner 
zone  electron  data  were  obtained.  Figure  8  shows  the  equatorial  flux  contours  as  a  func¬ 
tion  of  L  and  energy  provided  by  AP8MIN.  Since  the  MIN  model  predicts  slightly  more 
flux  than  the  MAX  model,  it  can  be  used  during  solar  maximum  or  for  long  term  mis¬ 
sions  as  a  conservative  model.  One  major  problem  with  the  energetic  proton  models  is 
the  fact  that  they  are  organized  in  terms  of  B,L.  The  secular  variation  in  the  earth’s 


magnetic  field  (the  dipole  term  is  diminishing)  causes  the  energetic  proton  ensemble, 
which  is  nominally  very  stable,  to  be  carried  to  lower  altitude.  The  model  does  not  take 
into  account  the  increased  atmospheric  density  the  protons  will  encounter  at  the  lower 
altitude.  As  a  result,  if  calculations  are  made  with  the  magnetic  field  projected  well  into 
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Figure  8.  AP8MIN  Equatorial  Proton  Flux  Contours  [15]. 
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the  future  (more  than  ten  to  fifteen  years),  totally  invalid  results  are  obtained  for  low 
altitudes  such  as  the  Space  Station  orbit  [16].  Since  the  source  of  the  energetic  protons  is 
the  decay  of  energetic  neutrons  produced  in  the  upper  atmosphere  by  cosmic  rays,  the 
geometry  of  the  production  process  relative  to  the  atmosphere  will  not  change.  The 
future  configuration  of  the  inner  zone  proton  belt  probably  will  not  change  relative  to  the 
present  configuration,  provided  both  are  described  in  terms  of  L  and  B/B0,  where  B  is  the 
magnetic  field  at  the  point  in  question  and  B0  is  the  equatorial  intensity  on  the  same  field 
tine.  One  will  almost  certainly  get  a  more  accurate  result  for  a  calculation  of  the  proton 
environment  for  Space  Station  in  the  year  2025  by  making  the  calculation  with  the 
present  field  model  than  by  extrapolating  the  field  35  years  into  the  future. 

Solar  Flare  Proton  Models 

Solar  flare  protons  are  treated  separately  because  they  are  transient  and  show  up 
in  the  outer  zone  only  in  relatively  small  numbers.  Their  major  impact  is  in  the  case  of 
low  altitude  polar-orbiting  satellites  which  normally  receive  relatively  tittle  particle  irra- 
diation-and  almost  all  of  that  in  short  infrequent  traversals  of  the  South  Atlantic 
Anomaly  region.  For  some  satellites,  the  major  portion  of  the  particle  environment  they 
encounter  during  their  lifetime  may  occur  during  one  or  two  major  solar  flare  particle 
events,  each  lasting  only  a  few  days.  In  a  typical  solar  cycle  of  1 1  years,  90%  of  the  en¬ 
ergetic  proton  fluence  is  the  result  of  a  single  anomalously  large  flare.  The  integrated  flu- 
ence  over  the  polar  caps  from  the  one  event  can  be  of  the  order  of  2  x  lO'Vcmz  P+  >  30 
MeV.  Averaged  over  a  solar  cycle,  the  annual  fluence  is  about  5  x  lOVcntf  p+  >  30  MeV 
over  the  polar  cap.  At  100  MeV,  these  numbers  are  about  a  factor  of  30  lower.  A  low 
altitude  polar  orbiting  satellite  spends  roughly  40%  of  its  time  at  latitudes  where  these 
protons  can  gain  access.  Thus  such  a  satellite  may  receive  virtually  all  of  its  energetic 
proton  dose  from  solar  flares  rather  than  from  the  magnetospherically-trapped  protons. 
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FUTURE  MODELLING  REQUIREMENTS 


The  major  efforts  in  modelling  must  be  directed  toward  the  electron  environ¬ 
ments,  since  the  energetic  inner  zone  protons  are  relatively  stable  and  well  understood 
and  the  outer  zone  protons  can  be  described  well  statistically.  Within  limits,  their  evolu¬ 
tion  after  a  storm  is  also  generally  predictable.  This  is  primarily  due  to  the  relative  ease 
with  which  proton  measurements  can  be  made.  Also,  the  studies  of  the  physics  of  the 
protons  have  been  more  productive.  The  CRRES  mission  [12]  is  expected  to  provide  ad¬ 
ditional  high-quality  data  for  the  proton  dynamics  modelling. 

While  a  large  amount  of  data  is  available  which  has  not  yet  been  used  to  update 
the  various  trapped  radiation  models,  there  are  significant  gaps  in  the  electron  data  bases. 
A  major  reason  for  this  is  the  difficulty  of  separating  the  signal  from  a  small  flux  of  these 
particles  from  the  large  background  of  protons  and  cosmic  rays.  One  of  these  gaps  was 
mentioned  previously:  No  useful  data  base  has  ever  been  obtained  in  the  outer  zone  for 
electrons  above  300  keV  between  about  L=2.8  and  geosynchronous  orbit  at  geomagnetic 
latitudes  below  about  20’.  CRRES  will  lower  this  limit  to  10’.  Above  2.8  MeV  no  data 
are  available  for  altitudes  above  8000  km  except  for  some  geosynchronous  orbit  mea¬ 
surements  [14].  In  the  inner  zone  above  2  MeV,  no  usable  measurements  are  available 
because  of  tremendous  penetrating  proton  background  problems  in  detectors.  CRRES 
will  not  be  able  to  furnish  this  type  of  data,  either,  because  its  primary  radiation  mission 
is  an  engineering  one  and  electrons  above  2  MeV  in  the  inner  zone  can  be  ignored.  (If 
you  have  problems  measuring  them,  they  are  unlikely  to  cause  a  problem  on  operational 
vehicles.) 

One  major  area  which  has  not  been  modelled  is  the  dynamics  of  the  radiation 
belts  in  response  to  magnetic  activity.  One  report  [17]  attempts  to  correlate  the  response 
of  the  outer  zone  to  the  magnetic  index  Kp,  but  the  particle  data  were  all  obtained  at 
geosynchronous  orbit  and  the  field  lines  that  are  represented  in  the  Kp  index  are  high  lat¬ 
itude  field  lines  which  thread  through  the  geosynchronous  region.  Hence,  that  study  is 
limited  in  validity  to  the  geosynchronous  region.  Low  altitude  measurements  indicate 
that  in  general  the  outer  zone  does  not  correlate  with  Kp,  except  incidentally  when  major 


magnetic  storms  occur,  but  does  have  some  relationship  with  Dsl,  the  low  latitude  index 
mentioned  previously.  The  intention  is  to  get  the  type  of  data  base  from  CRRES  that  one 
needs  to  address  the  dynamics  modelling. 

In  order  to  be  useful,  the  dynamics  models  have  to  address  three  issues  quantita¬ 
tively :  a)  The  prediction  of  magnetic  storms  which  produce  increases  in  the  outer  zone 
energetic  electron  and  proton  populations;  b)  The  energy  spectra  and  L  profiles  of  the 
fluxes  as  a  function  of  the  magnetic  storm  parameters;  and  c)  The  evolution  of  the  dis¬ 
tribution  as  a  function  of  time  (and  magnetic  activity)  post-storm.  The  first  of  these  re¬ 
quires  a  better  understanding  of  the  solar-terrestrial  coupling  physics  and  may  have  to 
await  completion  of  the  Global  Geophysics  Program  (aka  International  Solar-Terrestrial 
Physics  program)  which  is  scheduled  to  start  launching  satellites  in  the  1993  time  period 
[12j.  A  first  step  for  b)  and  c)  would  be  to  make  a  statistical  model  of  particle  storm  and 
post-storm  behavior  from  a  large  data  base.  Again,  CRRES  is  designed  to  provide  the 
data  base  for  such  an  effort  and  current  plans  include  developing  such  a  statistical  model. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation’s  Laboratory  Operations  conducts 
experimental  and  theoretical  Investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  Investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  Including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -field-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements, 
diagnostics  and  radlometry,  mlcrowave/mllllmeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena ,  space  communication  systems . 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
Instrumentation. 


